Molecular mechanisms of self-incompatibility in Brassicaceae and Solanaceae.

Self-incompatibility (SI) is a mechanism for preventing self-fertilization in flowering plants. SI is controlled by a single S-locus with multiple haplotypes (S-haplotypes). When the pistil and pollen share the same S-haplotype, the pollen is recognized as self and rejected by the pistil. This review introduces our research on Brassicaceae and Solanaceae SI systems to identify the S-determinants encoded at the S-locus and uncover the mechanisms of self/nonself-discrimination and pollen rejection. The recognition mechanisms of SI systems differ between these families. A self-recognition system is adopted by Brassicaceae, whereas a collaborative nonself-recognition system is used by Solanaceae. Work by our group and subsequent studies indicate that plants have evolved diverse SI systems.

The phylogeographic journey of a plant species from lowland to highlands during the Pleistocene.

Phylogeographic history refers to how species evolve and diversify in response to historical, ecological, and demographic factors. The climate fluctuation during the Pleistocene period marked a crucial time in shaping many species' distribution and genetic structure, particularly those from southern South American grasslands. This work investigated the phylogeographic history of a highland grassland, Petunia altiplana T. Ando & Hashim. (Solanaceae), its diversity, and geographic distribution using a population genomic approach based on RAD-seq data. Our results indicated that, during the Pleistocene, when the grasslands expanded to highlands, the lowland populations of P. altiplana reached the higher open fields, enlarging their geographic distribution. We found that the P. altiplana genetic diversity followed the geographic division into eastern (E) and western (WE) population groups, with a subtle division in the E group regarding the Pelotas River headwater. The results also showed that isolation by distance was the main divergence pattern, with elevation playing a pivotal role in shaping WE and E groups. Our findings indicated that lowland-adapted populations quickly colonized highlands during the late Pleistocene.

Chromosome level genome assembly of endangered medicinal plant Anisodus tanguticus.

Anisodus tanguticus is a medicinal herb that belongs to the Anisodus genus of the Solanaceae family. This endangered herb is mainly distributed in Qinghai-Tibet Plateau. In this study, we combined the Illumina short-read, Nanopore long-read and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to de novo assemble the A. tanguticus genome. A high-quality chromosomal-level genome assembly was obtained with a genome size of 1.26 Gb and a contig N50 of 25.07 Mb. Of the draft genome sequences, 97.47% were anchored to 24 pseudochromosomes with a scaffold N50 of 51.28 Mb. In addition, 842.14 Mb of transposable elements occupying 66.70% of the genome assembly were identified and 44,252 protein-coding genes were predicted. The genome assembly of A. tanguticus will provide genetic repertoire to understand the adaptation strategy of Anisodus species in the plateau, which will further promote the conservation of endangered A. tanguticus resources.

A chromosome-level genome assembly of anesthetic drug-producing Anisodus acutangulus provides insights into its evolution and the biosynthesis of tropane alkaloids.

Tropane alkaloids (TAs), which are anticholinergic agents, are an essential class of natural compounds, and there is a growing demand for TAs with anesthetic, analgesic, and spasmolytic effects. Anisodus acutangulus (Solanaceae) is a TA-producing plant that was used as an anesthetic in ancient China. In this study, we assembled a high-quality, chromosome-scale genome of A. acutangulus with a contig N50 of 7.4 Mb. A recent whole-genome duplication occurred in A. acutangulus after its divergence from other Solanaceae species, which resulted in the duplication of ADC1 and UGT genes involved in TA biosynthesis. The catalytic activities of H6H enzymes were determined for three Solanaceae plants. On the basis of evolution and co-expressed genes, AaWRKY11 was selected for further analyses, which revealed that its encoded transcription factor promotes TA biosynthesis by activating AaH6H1 expression. These findings provide useful insights into genome evolution related to TA biosynthesis and have potential implications for genetic manipulation of TA-producing plants.

An evolutionary view of vanillylamine synthase pAMT, a key enzyme of capsaicinoid biosynthesis pathway in chili pepper.

Pungent capsaicinoid is synthesized only in chili pepper (Capsicum spp.). The production of vanillylamine from vanillin is a unique reaction in the capsaicinoid biosynthesis pathway. Although putative aminotransferase (pAMT) has been isolated as the vanillylamine synthase gene, it is unclear how Capsicum acquired pAMT. Here, we present a phylogenetic overview of pAMT and its homologs. The Capsicum genome contained 5 homologs, including pAMT, CaGABA-T1, CaGABA-T3, and two pseudogenes. Phylogenetic analysis indicated that pAMT is a member of the Solanaceae cytoplasmic GABA-Ts. Comparative genome analysis found that multiple copies of GABA-T exist in a specific Solanaceae genomic region, and the cytoplasmic GABA-Ts other than pAMT are located in the region. The cytoplasmic GABA-T was phylogenetically close to pseudo-GABA-T harboring a plastid transit peptide (pseudo-GABA-T3). This suggested that Solanaceae cytoplasmic GABA-Ts occurred via duplication of a chloroplastic GABA-T ancestor and subsequent loss of the plastid transit signal. The cytoplasmic GABA-T may have been translocated from the specific Solanaceae genomic region during Capsicum divergence, resulting in the current pAMT locus. A recombinant protein assay demonstrated that pAMT had higher vanillylamine synthase activity than those of other plant GABA-Ts. pAMT was expressed exclusively in the placental septum of mature green fruit, whereas tomato orthologs SlGABA-T2/4 exhibit a ubiquitous expression pattern in plants. These findings suggested that both the increased catalytic efficiency and transcriptional changes in pAMT may have contributed to establish vanillylamine synthesis in the capsaicinoid biosynthesis pathway. This study provides insights into the establishment of pungency in the evolution of chili peppers.

Multiple independent losses of the biosynthetic pathway for two tropane alkaloids in the Solanaceae family.

Hyoscyamine and scopolamine (HS), two valuable tropane alkaloids of significant medicinal importance, are found in multiple distantly related lineages within the Solanaceae family. Here we sequence the genomes of three representative species that produce HS from these lineages, and one species that does not produce HS. Our analysis reveals a shared biosynthetic pathway responsible for HS production in the three HS-producing species. We observe a high level of gene collinearity related to HS synthesis across the family in both types of species. By introducing gain-of-function and loss-of-function mutations at key sites, we confirm the reduced/lost or re-activated functions of critical genes involved in HS synthesis in both types of species, respectively. These findings indicate independent and repeated losses of the HS biosynthesis pathway since its origin in the ancestral lineage. Our results hold promise for potential future applications in the artificial engineering of HS biosynthesis in Solanaceae crops.

Genome-Wide Analysis of MBF1 Family Genes in Five Solanaceous Plants and Functional Analysis of SlER24 in Salt Stress.

Multiprotein bridging factor 1 (MBF1) is an ancient family of transcription coactivators that play a crucial role in the response of plants to abiotic stress. In this study, we analyzed the genomic data of five Solanaceae plants and identified a total of 21 MBF1 genes. The expansion of MBF1a and MBF1b subfamilies was attributed to whole-genome duplication (WGD), and the expansion of the MBF1c subfamily occurred through transposed duplication (TRD). Collinearity analysis within Solanaceae species revealed collinearity between members of the MBF1a and MBF1b subfamilies, whereas the MBF1c subfamily showed relative independence. The gene expression of SlER24 was induced by sodium chloride (NaCl), polyethylene glycol (PEG), ABA (abscisic acid), and ethrel treatments, with the highest expression observed under NaCl treatment. The overexpression of SlER24 significantly enhanced the salt tolerance of tomato, and the functional deficiency of SlER24 decreased the tolerance of tomato to salt stress. SlER24 enhanced antioxidant enzyme activity to reduce the accumulation of reactive oxygen species (ROS) and alleviated plasma membrane damage under salt stress. SlER24 upregulated the expression levels of salt stress-related genes to enhance salt tolerance in tomato. In conclusion, this study provides basic information for the study of the MBF1 family of Solanaceae under abiotic stress, as well as a reference for the study of other plants.

Flower Development in the Solanaceae.

Flower development is the process leading from a reproductive meristem to a mature flower with fully developed floral organs. This multi-step process is complex and involves thousands of genes in intertwined regulatory pathways; navigating through the FLOR-ID website will give an impression of this complexity and of the astonishing amount of work that has been carried on the topic (Bouché et al., Nucleic Acids Res 44:D1167-D1171, 2016). Our understanding of flower development mostly comes from the model species Arabidopsis thaliana, but numerous other studies outside of Brassicaceae have helped apprehend the conservation of these mechanisms in a large evolutionary context (Moyroud and Glover, Curr Biol 27:R941-R951, 2017; Smyth, New Phytol 220:70-86, 2018; Soltis et al., Ann Bot 100:155-163, 2007). Integrating additional species and families to the research on this topic can only advance our understanding of flower development and its evolution.In this chapter, we review the contribution that the Solanaceae family has made to the comprehension of flower development. While many of the general features of flower development (i.e., the key molecular players involved in flower meristem identity, inflorescence architecture or floral organ development) are similar to Arabidopsis, our main objective in this chapter is to highlight the points of divergence and emphasize specificities of the Solanaceae. We will not discuss the large topics of flowering time regulation, inflorescence architecture and fruit development, and we will restrict ourselves to the mechanisms included in a time window after the floral transition and before the fertilization. Moreover, this review will not be exhaustive of the large amount of work carried on the topic, and the choices that we made to describe in large details some stories from the literature are based on the soundness of the functional work performed, and surely as well on our own preferences and expertise.First, we will give a brief overview of the Solanaceae family and some of its specificities. Then, our focus will be on the molecular mechanisms controlling floral organ identity, for which extended functional work in petunia led to substantial revisions to the famous ABC model. Finally, after reviewing some studies on floral organ initiation and growth, we will discuss floral organ maturation, using the examples of the inflated calyx of the Chinese lantern Physalis and petunia petal pigmentation.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes.

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

Evolution and Expression of the Meprin and TRAF Homology Domain-Containing Gene Family in Solanaceae.

Meprin and TRAF homology (MATH)-domain-containing proteins are pivotal in modulating plant development and environmental stress responses. To date, members of the MATH gene family have been identified only in a few plant species, including Arabidopsis thaliana, Brassica rapa, maize, and rice, and the functions of this gene family in other economically important crops, especially the Solanaceae family, remain unclear. The present study identified and analyzed 58 MATH genes from three Solanaceae species, including tomato (Solanum lycopersicum), potato (Solanum tuberosum), and pepper (Capsicum annuum). Phylogenetic analysis and domain organization classified these MATH genes into four groups, consistent with those based on motif organization and gene structure. Synteny analysis found that segmental and tandem duplication might have contributed to MATH gene expansion in the tomato and the potato, respectively. Collinearity analysis revealed high conservation among Solanaceae MATH genes. Further cis-regulatory element prediction and gene expression analysis showed that Solanaceae MATH genes play essential roles during development and stress response. These findings provide a theoretical basis for other functional studies on Solanaceae MATH genes.

Genomic and structural basis for evolution of tropane alkaloid biosynthesis.

The tropane alkaloids (TAs) cocaine and hyoscyamine have been used medicinally for thousands of years. To understand the evolutionary origins and trajectories of serial biosynthetic enzymes of TAs and especially the characteristic tropane skeletons, we generated the chromosome-level genome assemblies of cocaine-producing Erythroxylum novogranatense (Erythroxylaceae, rosids clade) and hyoscyamine-producing Anisodus acutangulus (Solanaceae, asterids clade). Comparative genomic and phylogenetic analysis suggested that the lack of spermidine synthase/N-methyltransferase (EnSPMT1) in ancestral asterids species contributed to the divergence of polyamine (spermidine or putrescine) methylation in cocaine and hyoscyamine biosynthesis. Molecular docking analysis and key site mutation experiments suggested that ecgonone synthases CYP81AN15 and CYP82M3 adopt different active-site architectures to biosynthesize the same product ecgonone from the same substrate in Erythroxylaceae and Solanaceae. Further synteny analysis showed different evolutionary origins and trajectories of CYP81AN15 and CYP82M3, particularly the emergence of CYP81AN15 through the neofunctionalization of ancient tandem duplication genes. The combination of structural biology and comparative genomic analysis revealed that ecgonone methyltransferase, which is responsible for the biosynthesis of characteristic 2-substituted carboxymethyl group in cocaine, evolved from the tandem copies of salicylic acid methyltransferase by the mutations of critical E216 and S153 residues. Overall, we provided strong evidence for the independent origins of serial TA biosynthetic enzymes on the genomic and structural level, underlying the chemotypic convergence of TAs in phylogenetically distant species.

Anatomical and biophysical basis for graft incompatibility within the Solanaceae.

Interspecies grafting is an economically relevant technique that allows beneficial shoot and root combinations from separate species to be combined. One hypothesis for the basis of graft compatibility revolves around taxonomic relatedness. To test how phylogenetic distance affects interspecific graft compatibility within the economically important Solanaceae subfamily, Solanoideae, we characterized the anatomical and biophysical integrity of graft junctions between four species: tomato (Solanum lycopersicum), eggplant (Solanum melongena), pepper (Capsicum annuum), and groundcherry (Physalis pubescens). We analyzed the survival, growth, integrity, and cellular composition of the graft junctions. Utilizing various techniques, we were able to quantitatively assess compatibility among the interspecific grafts. Even though most of our graft combinations could survive, we show that only intrageneric combinations between tomato and eggplant are compatible. Unlike incompatible grafts, the formation of substantial vascular reconnections between tomato and eggplant in the intrageneric heterografts likely contributed to biophysically stable grafts. Furthermore, we identified 10 graft combinations that show delayed incompatibility, providing a useful system to pursue deeper work into graft compatibility. This work provides new evidence that graft compatibility may be limited to intrageneric combinations within the Solanoideae subfamily. Further research amongst additional Solanaceous species can be used to test the extent to which our hypothesis applies to this family.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants.

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family.

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

The SlDOG1 Affect Biosynthesis of Steroidal Glycoalkaloids by Regulating GAME Expression in Tomato.

Steroidal alkaloids (SAs) and steroidal glycoalkaloids (SGAs) are common constituents of plant species belonging to the Solanaceae family. However, the molecular mechanism regulating the formation of SAs and SGAs remains unknown. Here, genome-wide association mapping was used to elucidate SA and SGA regulation in tomatoes: a SlGAME5-like glycosyltransferase (Solyc10g085240) and the transcription factor SlDOG1 (Solyc10g085210) were significantly associated with steroidal alkaloid composition. In this study, it was found that rSlGAME5-like can catalyze a variety of substrates for glycosidation and can catalyze SA and flavonol pathways to form O-glucoside and O-galactoside in vitro. The overexpression of SlGAME5-like promoted α-tomatine, hydroxytomatine, and flavonol glycoside accumulation in tomatoes. Furthermore, assessments of natural variation combined with functional analyses identified SlDOG1 as a major determinant of tomato SGA content, which also promoted SA and SGA accumulation via the regulation of GAME gene expression. This study provides new insights into the regulatory mechanisms underlying SGA production in tomatoes.

Establishing Physalis as a Solanaceae model system enables genetic reevaluation of the inflated calyx syndrome.

The highly diverse Solanaceae family contains several widely studied models and crop species. Fully exploring, appreciating, and exploiting this diversity requires additional model systems. Particularly promising are orphan fruit crops in the genus Physalis, which occupy a key evolutionary position in the Solanaceae and capture understudied variation in traits such as inflorescence complexity, fruit ripening and metabolites, disease and insect resistance, self-compatibility, and most notable, the striking inflated calyx syndrome (ICS), an evolutionary novelty found across angiosperms where sepals grow exceptionally large to encapsulate fruits in a protective husk. We recently developed transformation and genome editing in Physalis grisea (groundcherry). However, to systematically explore and unlock the potential of this and related Physalis as genetic systems, high-quality genome assemblies are needed. Here, we present chromosome-scale references for P. grisea and its close relative Physalis pruinosa and use these resources to study natural and engineered variations in floral traits. We first rapidly identified a natural structural variant in a bHLH gene that causes petal color variation. Further, and against expectations, we found that CRISPR-Cas9-targeted mutagenesis of 11 MADS-box genes, including purported essential regulators of ICS, had no effect on inflation. In a forward genetics screen, we identified huskless, which lacks ICS due to mutation of an AP2-like gene that causes sepals and petals to merge into a single whorl of mixed identity. These resources and findings elevate Physalis to a new Solanaceae model system and establish a paradigm in the search for factors driving ICS.

Down, then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species, Nicotiana section Suaveolentes (Solanaceae).

BACKGROUND AND AIMS: The extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia. METHODS: We analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework. KEY RESULTS: RADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15-18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids. CONCLUSIONS: The results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.

Phylogenetic Analyses of Some Key Genes Provide Information on Pollinator Attraction in Solanaceae.

Floral syndromes are known by the conserved morphological traits in flowers associated with pollinator attraction, such as corolla shape and color, aroma emission and composition, and rewards, especially the nectar volume and sugar concentration. Here, we employed a phylogenetic approach to investigate sequences of genes enrolled in the biosynthetic pathways responsible for some phenotypes that are attractive to pollinators in Solanaceae genomes. We included genes involved in visible color, UV-light response, scent emission, and nectar production to test the hypothesis that these essential genes have evolved by convergence under pollinator selection. Our results refuted this hypothesis as all four studied genes recovered the species' phylogenetic relationships, even though some sites were positively selected. We found differences in protein motifs among genera in Solanaceae that were not necessarily associated with the same floral syndrome. Although it has had a crucial role in plant diversification, the plant-pollinator interaction is complex and still needs further investigation, with genes evolving not only under the influence of pollinators, but by the sum of several evolutionary forces along the speciation process in Solanaceae.

Complete Plastome of Physalis angulata var. villosa, Gene Organization, Comparative Genomics and Phylogenetic Relationships among Solanaceae.

Physalis angulata var. villosa, rich in withanolides, has been used as a traditional Chinese medicine for many years. To date, few extensive molecular studies of this plant have been conducted. In the present study, the plastome of P. angulata var. villosa was sequenced, characterized and compared with that of other Physalis species, and a phylogenetic analysis was conducted in the family Solanaceae. The plastome of P. angulata var. villosa was 156,898 bp in length with a GC content of 37.52%, and exhibited a quadripartite structure typical of land plants, consisting of a large single-copy (LSC, 87,108 bp) region, a small single-copy (SSC, 18,462 bp) region and a pair of inverted repeats (IR: IRA and IRB, 25,664 bp each). The plastome contained 131 genes, of which 114 were unique and 17 were duplicated in IR regions. The genome consisted of 85 protein-coding genes, eight rRNA genes and 38 tRNA genes. A total of 38 long, repeat sequences of three types were identified in the plastome, of which forward repeats had the highest frequency. Simple sequence repeats (SSRs) analysis revealed a total of 57 SSRs, of which the T mononucleotide constituted the majority, with most of SSRs being located in the intergenic spacer regions. Comparative genomic analysis among nine Physalis species revealed that the single-copy regions were less conserved than the pair of inverted repeats, with most of the variation being found in the intergenic spacer regions rather than in the coding regions. Phylogenetic analysis indicated a close relationship between Physalis and Withania. In addition, Iochroma, Dunalia, Saracha and Eriolarynx were paraphyletic, and clustered together in the phylogenetic tree. Our study published the first sequence and assembly of the plastome of P. angulata var. villosa, reported its basic resources for evolutionary studies and provided an important tool for evaluating the phylogenetic relationship within the family Solanaceae.

Elucidation of tropane alkaloid biosynthesis in Erythroxylum coca using a microbial pathway discovery platform.

Tropane alkaloids (TAs) are heterocyclic nitrogenous metabolites found across seven orders of angiosperms, including Malpighiales (Erythroxylaceae) and Solanales (Solanaceae). Despite the well-established euphorigenic properties of Erythroxylaceae TAs like cocaine, their biosynthetic pathway remains incomplete. Using yeast as a screening platform, we identified and characterized the missing steps of TA biosynthesis in Erythroxylum coca. We first characterize putative E. coca polyamine synthase- and amine oxidase-like enzymes in vitro, in yeast, and in planta to show that the first tropane ring closure in Erythroxylaceae occurs via bifunctional spermidine synthase/N-methyltransferases and both flavin- and copper-dependent amine oxidases. We next identify a SABATH family methyltransferase responsible for the 2-carbomethoxy moiety characteristic of Erythroxylaceae TAs and demonstrate that its coexpression with methylecgonone reductase in yeast engineered to express the Solanaceae TA pathway enables the production of a hybrid TA with structural features of both lineages. Finally, we use clustering analysis of Erythroxylum transcriptome datasets to discover a cytochrome P450 of the CYP81A family responsible for the second tropane ring closure in Erythroxylaceae, and demonstrate the function of the core coca TA pathway in vivo via reconstruction and de novo biosynthesis of methylecgonine in yeast. Collectively, our results provide strong evidence that TA biosynthesis in Erythroxylaceae and Solanaceae is polyphyletic and that independent recruitment of unique biosynthetic mechanisms and enzyme classes occurred at nearly every step in the evolution of this pathway.